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Abstract 



We investigate the role played by subcritical bubbles at the onset of the 
electroweak phase transition. Treating the configuration modelling the ther- 
mal fluctuations around the homogeneous zero configuration of the Higgs 
field as a stochastic variable, we describe its dynamics by a phenomenolog- 
ical Langevin equation. This approach allows to properly take into account 
both the effects of the thermal bath on the system: a systematic dyssi- 
pative force, which tends to erase out any initial subcritical configuration, 
and a random stochastic force responsible for the fluctuations. We show 
that the contribution to the variance ((<p 2 {t)))v i n a given volume V from 
any initial subcritical configuration is quickly damped away and that, in 
the limit of long times, (((fi 2 (t)))v approaches its equilibrium value provided 
by the stochastic force and independent from the viscosity coefficient, as 
predicted by the fluctuation-dissipation theorem. In agreement with some 
recent claims, we conclude that thermal fluctuations do not affect the nu- 
cleation of critical bubbles at the onset of the electroweak phase transition 
making electroweak baryogenesis scenarios still a viable possibility to explain 
the primordial baryon asymmetry in the Universe. 
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1. Introduction 

Nucleation of critical bubbles during a first order electroweak phase tran- 
sition has received much attention since the discovery of the possibility for 
electroweak baryogenesis [0] . Indeed, one of the basic ingredients for the gen- 
eration of the baryon asymmetry (apart from the requirement of baryon- and 
CP-violating interactions) is the presence of an out-of-equilibrium state || 
which, during the first order electroweak phase transition with supercooling, 
is attained by critical bubbles expanding in the thermal bath of the unbroken 
phase. 

However, less attention has been paid to the environment where the crit- 
ical bubble nucleation is supposed to occur. Since critical bubbles have a 
finite size, phase transitions are highly local phenomena. Fluctuations of 
the Higgs scalar field <f> around the origin = with spatial correlations 
comparable to the critical bubble size may be expected to be important for 
bubble nucleation. Also, if thermal fluctuations are too large, any perturba- 
tive scheme could break down. In such a situation, the prediction of a first 
order phase transition becomes suspect, with the possibility that the entire 
scenario of electroweak baryogenesis might be invalidated. 

On the other hand, while the presence of thermal fluctuations in any 
thermodynamical system is undisputed, their role in the dynamics of weakly 
first order phase transitions is still controversial. 

In ref. 0, it was first conjectured that statistical fluctuations around 
equilibrium may be described by spherically symmetric configurations roughly 
extending over a correlation volume, where the correlation length is given by 
the inverse temperature dependent mass of the Higgs field: £(T) = m _1 (T). 

These fluctuations are referred to as subcritical bubbles. Their amplitude 
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was estimated in ref. 0, where it was concluded that they are dominant if 
the Higgs mass M# is larger than ~ 80 GeV, whence the fraction which the 
asymmetric vacuum occupies in the neighborhood of the critical tempera- 
ture becomes of the order of unity. Therefore it was concluded that critical 
bubbles cannot be generated due to the inhomogeneities of the background 
field. 

In ref. ||, however, the continuous disappearance of the subcritical bub- 
bles was not accounted for. This can happen in two ways: the subcritical 
bubbles, being unstable configurations, tend to shrink; they are also sub- 
ject to constant thermal bombardment so that they may disappear simply 
because of thermal noise. 

Following general principles and estimating the scalar field two-point 
function, i.e. the variance of a gaussian fluctuation distribution, computed 
in a correlation volume, Dine et al. [|J and Anderson [||] have later argued 
that subcritical bubbles do not affect the nucleation of critical bubbles in an 
appreciably way. 

Subsequently, Gelmini and Gleiser rekindled the issue adopting a dif- 
ferent point of view based on modelling by a set of Boltzmann equations the 
evolution with time t of the number density n(R,t) of subcritical bubbles 
with a certain radius R. 

Under a specific assumption about the form of the destruction rate due 
to thermal noise, they found that for Higgs masses below ~ 55 GeV the 
approach to equilibrium is dominated by shrinking. Unfortunately, for the 
interesting range of Higgs masses dictated by the experimental constraints 
coming from LEP, Mh > 60 GeV 0, their analysis is inconclusive since the 
approximations adopted break down. 

Recently, Enqvist et al. |5J re-estimated the amplitude, average size and 
formation rate of subcritical bubbles taking into account the crucial role 
played by thermalization. Their starting point was the the observation that 
at the microscopic level, the true origin of the dominant thermal fluctuation is 
the perpetual creation and annihilation of spherical subcritical bubbles. Thus 
one should identify the typical amplitude and size of these bubbles with the 
values estimated by a statistical ensemble averaging, instead of assuming, for 
instance, the size fixed and equal to the correlation length £(T). 
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They also observed that a large subcritical bubble should resemble the 
critical one around the critical temperature: when R increases, the form 
of the subcritical bubble should deform smoothly to reproduce the critical 
configuration when R = R c , R c being the critical radius. Thus, it was found 
that the average size (R) of the subcritical bubbles is much larger than the 
correlation length, and that the average amplitude {(fo) 1 ^ 2 at the core is much 
smaller than any previous estimate. 

These results led the authors of ref. || to conclude that thermal fluctu- 
ations do not hinder the electroweak first order phase transition. Of course, 
treating the size and the amplitude of subcritical bubbles as statistical de- 
grees of freedom of an equilibrium ensamble at a fixed temperature T does not 
allow to follow the time evolution of such degrees of freedom. Thus, since in 
the approach of ref. || there is no dynamics involved, (R) and {4>o) 1 ^ 2 should 
be interpretated as the most probable initial conditions for any configuration 
describing a subcritical bubble generated at the time t = 0. 

Analogous conclusions to those of ref. || have been very recently obtained 
by Bettencourt M, who, along the same line followed by Hindmarsh and 
Rivers [ IT3 1 for the A0 4 theory, computed the probability for fluctuations of 
the Standard Model Higgs field, averaged over a given spatial scale, to exceed 
a specified value. He found that the probability for the Higgs scalar field to 
fluctuate from the symmetric to the asymmetric minimum before the latter 
becomes stable is very small for Higgs masses of order of those of the 
and Z° bosons, whereas the converse is more likely. 

The aim of the present paper is to take a further step in the investigation 
of the role played by subcritical bubbles during the onset of the electroweak 
phase transition. In particular, we wish to present an analytic treatment of 
the dynamics of subcritical bubbles, described by a suitably coarse-grained 
configuration (x, t), in a random environment. 

Indeed, thermal fluctuations are the manifestation of the interaction be- 
tween the system and the surrounding environment. The peculiarity of the 
thermal bath is that a virtually infinite number of degrees of freedom takes 
part in the exchanges of energy with the system. A familiar example |TTJ are 
the molecules of a gas or a fluid in which the system (the brownian particle) 
is embedded. It is well known that the impacts between the the system and 
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the surrounding particles cause two effects: a systematic dissipative force 
(friction) and a random force, responsible for the fluctuations. 

To describe the dynamics of subcritical fluctuations in a thermal bath 
around the homogeneous configuration = 0, we will then assume a classical 
Markovian Langevin equation 

0(x, t) - V 2 (x, t) + #(x, t) = - 9V{ g' T) + 7 (x, t) , (1) 

where V{<f>, T) is the potential associated to the the scalar field <fi (x, t) and 
1] is the viscosity coefficient which takes into account the dissipative effects. 

The fluctuations are modelled by introducing a stochastic additive force 
(noise term) 7 (x, t) characterized by a gaussian distribution with 

(7 CM)) = 0, 

( 7 (x,t) 7 (x',t')> = 2D5 3 (x-x')5(t-i'), (2) 

where D is the diffusion coefficient. 

The Langevin equation (1) describing the dynamics of a subcritical bubble 
in contact with a thermal bath is similar to the Kardar-Parisi-Zhang equation 
[ |T^j modelling the evolution of fluctuations around the average profile of a 
growing interface, which, in our case, is represented by the homogeneous 
configuration = 0. 

Since the dynamics is expected to bring the system into equilibrium with 
the thermal bath at long times, the Langevin equation describes correctly 
the long-time dynamics and reproduce the correct equilibrium behaviour if 
the fluctuation-dissipation condition D = rjT, expressing the common origin 
of dissipation and fluctuation, is satisfied [|TB|]. 

Equation (1) naturally incorporates different crucial informations: i) sub- 
critical bubbles are subject to constant thermal bombardment so that they 
can be rapidly thermalized and disappear with a relaxation rate rj (this was 
also accounted for in ref. || where, however, no dynamics was involved); 
ii) subcritical bubbles, being unstable configurations, tend to shrink; in) 
the crucial role played by the noise term 7 (x, t) responsible for fluctuations. 
It tends to contrast the damping term 7]<j) and determines the form of the 
fluctuation (0 2 (x, £)) at long times as well as its final equilibrium value. 

A few comments are in order here about our choice of eq. (1). In writing 
it we have assumed that the system is Markovian: the correlation time scale 
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for the noise is smaller than the typical relaxation time scale for the system. 
We have also assumed the noise to be additive. 

Recent works indicates that one should expect departures from the Langevin 
equation written above |H |, although the details are very sensitive to the 
model one starts with. For istance, in a A0 4 theory, the Langevin equation 
describing the dynamics of the long wavelength modes evolving in the bath 
formed by the short wavelength modes is characterized by a colored and 
multiplicative noise and by a space-time dependent viscosity coefficient (see 
Gleiser and Ramos in 0]). 



In this paper we will adopt eq. (1) as a first step. In fact, we do not 
expect that the nature of the thermal noise will change the final equilibrium 
properties of the system, even if the relaxation time-scales can be changed. 
Since the physical results will be related to the final equilibrium value of 
(0 2 (x, t)), we believe that they will be not affected by more complicated 
(even if more physical) representations of the coupling of the field </> (x, t) to 
the thermal bath. 

This point of view is also motivated by the independence of the statistical 
moments of the (x, t) field at long times from the viscosity coefficient r], as 
implicit in the fluctuation-dissipation theorem. 

For sake of clarity we have decided to address the issue of dynamics of 
subcritical bubbles in the electroweak theory after having described the same 
issue for simpler theories. 

The paper is then organized as follows. In Subsection 2.1 we describe in 
details, for a free scalar field theory in the absence of thermal environments, 
the dynamics of a configuration representing an initial departure from the 
minimum energy equilibrium configuration. This allows us to evaluate the 
typical shrinking lifetime r s h of such a configuration. 

In Subsection 2.2 we let the free scalar field be coupled to a thermal 
environment and follow the dynamics of the fluctuation at short and long 
times. We determine the asymptotic value of (0 2 (x, t)) in a given volume for 
t ^> 1/r], showing explicitly that it is independent of time and equal to its 
equilibrium value. 

In Section 3 we make use of the results obtained in Section 2 to describe, 
through a self-consistent Hartree approximation, the strength of the thermal 
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fluctuations for the A0 4 theory. 

In Section 4 we present our results for the electroweak theory, from which 
one can conclude that subcritical bubbles do not affect nucleations of critical 
bubbles in a appreciably way. 

We finally draw in Section 5 our conclusions and perspectives for future 
work. 



2. Free scalar field: linear dynamics 
2.1 Free evolution at zero temperature 

In this Section we want to analyze the evolution of a spherically-symmetric 
configuration in the linear regime and absence of coupling to the thermal 
bath. We already know from Derrick's theorem fl5f , which forbids in four 



dimensions the existence of static, finite-energy configurations for models 
containing a simple real scalar field, that such configurations must be un- 
stable and decay with a relaxation time r s h, that we now want to estimate 
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Given the quadratic mass term (m 2 > 0) 

VM = ^4>\ (3) 

the associated Klein-Gordon equation is given by 

0(x, f) - V 2 (x, t) + m 2 (x, t) =0. (4) 

We model the initial deviation at t = from the homogeneous minimum- 
energy configuration = by a gaussian shape of the type 

0(x,O) = 0o e-W 2 / 2 ^, 

0(x,O) = 0, (5) 

i.e. we imagine a deviation formed initially at rest with initial radius i?o 
and amplitude O at its core. Eq. (g) is most easily solved by examining the 
Fourier components 

(q, t) = J rf 3 x e iq - x (x, t) , (6) 
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which evolve according to 

0(q,i) + (|q| 2 + m 2 )0(q,i) = O. (7) 
The solution of equation (fil) is easily found to be 



(q, t) = A(q) cos (y 7 |q| 2 + m 2 t j + B(q) sin (y |q| 2 + m 2 tj . (8) 

The initial conditions described in eq. (||) fix -A(q) an d -B(q) 

A(q) = ( 27 r) 3 / 2 o ^e-l^^/ 2 ; 

B(q) = 0, (9) 

which, when plugged into eq. @ allow us to investigate the behaviour of the 
bubble's core with time 



(0,t) = ^% i?3 / o °°rf|q| |q| 2 e-l«/ 2 cos ( y^ 2 + m 2 tj 



(10) 



Since the integral is dominated by small values of |q|, |q| < \/2R 1 , we can 



approximate the argument of cos f y |q| + m 2 tj for Ro > \J~2vn 1 and write 
0(0, t) ~ \Jl<t>oRl J~d\q\ |q| 2 Re (e jm * e-l q l 2/? o/ 2 e ^lV2^ 



cos 



mt + | tan 1 (t/mRl 



= foRl l — —Hi L - (11) 



(■Ro + S) ' 

Thus, given an initial deviation from the stable configuration = with 
sufficiently large initial size Rq and in absence of coupling to a thermal bath, 
the amplitude at the core decays as t~ 3//2 and the unstable configuration 
shrinks by radiating away its initial energy after a lifetime ]16| 

r sh ~mRl>-. (12) 
m 



2.2 Coupling to the thermal bath 

The next step is to investigate how the dynamics of the spherically- 
symmetric fluctuations around the homogeneous configuration = changes 
when the system is coupled to a thermal bath. 
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The Langevin equation (1) reads 

0(x, t) - V 2 </> (x, t) + 7$(x, t) + m 2 4> (x, t) = 7 (x, t) . (13) 

One can imagine the above equation to emerge from an effective theory where 
the fields (for instance, scalars different from (x, t) or fermions) which 
(f> (x, t) is coupled to have been integrated out leaving an effective poten- 
tial approximated at high temperature by eq. (|3|). In such a case m 2 has 
to be understood as a function of the temperature, m 2 (T) ~ g 2 T 2 , g being 
the generic coupling constant of <f> (x, t) to the other fields of the underlying 
original theory. 

We wish to point out once more that the form of the Langevin equation 
arising from an effective theory is strongly model-dependent. Nevertheless, 
since we are interested in the equilibrium value of the fluctuation (0 2 (x, t)), 
we are confident that our result will not be affected by using eq. ([13]) to 
describe the dynamics of subcritical bubbles. 

Again, eq. ([TS]) is most easily solved by examining the Fourier components 
of the field (x, t), see eq. (|J). If we define 

7(q,t) = / d 3 xe^ 7 (x,f), (14) 

the Fourier transformed noise, the equation to solve is 

I (q, t) + 4 (q, t) + (|q| 2 + m 2 ) $ (q, t) = j{q, t). (15) 

We define the function \ as 

X (q) =77 2 -4(|q| 2 + m 2 ), (16) 

and assume 

7]<2m. (17) 

The case r] > 2 m, though straightforward, is computationally more involved. 
On the other hand, it is physically the least interesting case, since any oscil- 
latory behaviour is quickly dominated by the exponentially damped factors. 
Moreover, it obviously leads to the same results as far as the final equilibrium 
values are concerned. 

The solution of eq. ([15]) can be cast into the form 

0(q,t) =0 fto (q,t)+0 7 (q,t), (18) 
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where 



4>ho(ci,t) 



7 (q,t) 




Note that 7 (q, t) does not depend upon the initial conditions described by 
eq. (H), while they are embodied in 0/i O (q, t) through the amplitudes A(q) 
(given by eq. (|9])) and 

B(q) = ^=L=^(q). (20) 

vlx(q)l 

When computing the mean values (••) with respect to the gaussian measure 
of the noise 7(x, t), some care is needed. Indeed, as noted in the Introduction, 
one has to assume that the correlation time scale for the noise r 7 is smaller 
than the typical relaxation time of the free system r sh ~ m as determined 
in Subsect. 2.1. 

In the case under examination, rj < 2m, this is true only for configurations 
with very large initial radius Rq S> m , then, from now on we will assume 
this to be the caseQ, i.e. r 7 <C r s h- 

From eq. (|T9|) we see that fluctuations in each mode decay with a relax- 
ation time given by r 7 ~ 1/rj. This implies for the following decay law for 
the average values: 

(<Mq,t)> ~ exp(-77t/2), (21) 

as expected from the presence of the damping term in the Langevin equa- 
tion. This is similar to what happens for a particle thrown into a thermal 
bath with an initial velocity v®: after a sufficiently long time, and indepen- 
dently from the values of ■uo, the particle behaves like a brownian system 

3 Of course, in the opposite limit r\ > 2 m, the correlation time scale for the noise is 
much smaller than any other time scales and no additional condition on R is required. 
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with a configurational average (X) = and a nonvanishing fluctuation (X 2 ) 
determined by the surrounding noise. 

The two-point correlation function at equal time is given by 

(0(x,t)0(x',t)> = /^" ,q ' X /^*'^( qi t)fet)) 

= /lo (x,t)0 h , o (x / ,t) + (0(x,t)0(x / ,t)) 7 , (22) 

where we have defined 

<<Mx,t)<Kx',t)> 7 = / J^e*** f |^e^'- x '(0 7 (q,t)0 7 (q',t)). (23) 

Making use of the two-point correlation function for the Fourier transformed 
noise 

( 7 (q, t) 7 (q', f)> = 2 (2tt) 3 t? T <5 3 (q + q') *(f - t'), (24) 
eq. ( p^) becomes 

^ 3( 1 „-ia(x-xM 1 



(x, t )0(x', t )) 1 = 4T/^e-.M) |;((q) 



?7 2 + |x(q 



r/ 2 cos v /|x(q)|t-r/^|x(q)| sin^/|x(q)|t 

^ 2 + lx(q)l 

(25) 

We are now in the position to study the evolution with time of the magnitude 
of the fluctuation around the = state. A natural measure of it is given 
by the two-point connected Green's function, coarse-grained on a volume V, 
defined as 

«0 2 W»V = J d " X j rf 3 x'(0(x,t)0(x',t))/(x)/(x') 

= mt)))v + mtw v , (26) 

where the definition of ((4> 2 (t)))v° an d {{(ft 2 (t)))v can De easily read from eq. 
and 

'M = (27) 
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is a window function modelling the volume V = (4ir/3)R 3 over which we test 
the magnitude of the fluctuations^ 

The system as a whole is taken to have a volume much larger than V. 
If we denote by £ ~ 1/m the correlation length, we can parametrize the 
number of correlation volumes by writing R = (3£, where (3 is a positive 
number. Since, at any time, the field is correlated over some correlation 
length, it will be sufficient to take (3 > 1. 

From eq. (p6| ) we see that the variance {{<p 2 {t)))v receives two contri- 
butions: ((0 2 (^)))v r ° describing the contribution coming from the gaussian 
configuration formed at t = and {{4> 2 {t)))y expressing the contribution due 
to the noise force. As a consequence, the latter does not depend on the initial 
conditions. 

Let us now imagine that the gaussian configuration (|5]) represents a devi- 
ation from the equilibrium value formed at t = due to thermal fluctuations. 
To follow the dynamics of ((<p 2 {t)))v we can then identify two different time 
regimes: 

i) Short times: 

For t <C 1/f], none of the modes of ((<p 2 {t)))y° has yet relaxed, but {(<p 2 {t)))y° 
starts to decrease as t 2 , that is 



< 2 77#4t3 U + sis) ■ (28) 



Taking Rq ~ R — and /3 ^> 1, we see that the contribution to 

ho 
V ■ 



{{4> 2 (t)))v coming from the gaussian configuration formed at t = 0, ((4> 2 (t))) ho 



decreases as ~ t 2 4>lm 2 



4 The coefficient ^2/9ir in cq. (|27|) is chosen to make easier the comparison between 
our results in Section 4 and those in ref. || in the case of the fluctuations at the onset of 
the electroweak phase transition. 
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In the meanwhile, the contribution to ((<p 2 {t)))v from the noise force, 
{{4> 2 (t)))y, starts to increase as t 3 : 

2r]T a3 f d 3 q _i„i2 E?2 



1 VT- 



12 7T 3 / 2 R 3 



t s . (29) 



izj Long times: 



For t ^> I/77, all the modes have relaxed to their equilibrium values. The ini- 
tial gaussian configuration (|5]) has been rapidly thermalized with a relaxation 
time r 7 ~ I/77 so that {{4> 2 {t)))y' has quickly vanished. 

The situation is again similar to that for a brownian particle immersed 
in a thermal environment at t — with initial data X Q , v : the contribution 
to (X 2 (t)) due to the initial conditions is very rapidly dissipated away and 
taken over by the contribution from the noise. 

As a consequence, at long times {{(j) 2 (t)))y provides the only contribu- 
tion to ((<p 2 {t)))v which tends asymptotically towards the time-independent 
dynamical value 

((^ 2 )) d / n = / 7^-'^ 2 ^-. (30) 



• (2tt) 3 |q| 2 + m 2 " 

Note that the dependence from the the viscosity coefficient rj in the above 
expression is absent as implicit in the fluctuation-dissipation theorem: the 
system after a sufficiently long time relaxes towards its equilibrium state 
which is completely independent from 7]. 

It is readily verified that the above expression coincides with the equilib- 
rium value of the two-point connected Green's function, coarse-grained on a 
volume V, calculated in ref. [IC]. This is not surprising. We already know 



that dynamics is expected to bring the thermal fluctuations represented by 
the subcritical bubbles into equilibrium with the surrounding thermal bath 
at long times. 

The result expressed by eq. ([JU]) shows explicitly that the Langevin equa- 
tion correctly describes the long-time behaviour of ((<p 2 {t)))v reproducing its 
equilibrium value if the fluctuation-dissipation condition is satisfied, i. e 

lim^oo {{<P\t)))v = «0 2 » d / n = «0 2 »^ Q (31) 
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Note that eq. (^) still holds also in the case rj > 2 m. Namely, in the 
overdamped limit (large times or strong coupling to the environment), the 
description simplifies considerably, as the canonical momentum associated 
to <p (x, t) becomes a variable much faster than the field itself and can be 
considered fully thermalized in the adiabatic approximation. 

The above statement amounts to drop the inertia term 0(x, t) in eq. (|15|), 
leading to the equation 



V<f> (q, t) + (|qf + m 2 ) <P (q, t) = 7(q, t), (32) 

which reproduces the large time behaviour expressed by eq. (|30|). 

Assuming R = f3/m and noting that the gaussian function in eq. ([!(]) is 
substantially different from zero only for |q| 2 < m 2 //3 2 , we get 



\ - tan " Ik 



(33) 



The spreading of the field over a correlation volume = is then 

((0 2 »JT ^10- 2 mT. (34) 

For volumes V substantially larger than the correlation volume, ((4> 2 ))v rn — 
(m T/3 7r 2 (3 3 ) «mT. The coefficient 1/27T 2 plays a crucial role in reducing 
the field spreading, as already pointed out in refs. P, P|, |TT) . 



One could expect from these results that subcritical fluctuations will not 
play a significant role in the electroweak phase transition. Nevertheless, all 
the considerations made in this Section are limited to a linear system. In the 
next Section we wish to extend them, through a self-consistent Hartree ap- 
proximation, to a nonlinear system described by a A0 4 potential as a further 
step forwards the study of thermal fluctuations at the onset of electroweak 
phase transition. 



3. Self-interacting field: Hartree approximation 

The above description of the evolution of spherically-symmetric initial 
configurations in the presence of a thermal bath is based on the linearity of 
the Langevin equation. 
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Let us now consider the issue of thermal fluctuations around the homo- 
geneous configuration = for a nonlinear system with a potential 

V{</>) = + (35) 

where m 2 is assumed to be positive so that V((f>) posseses a unique minimum 
at 4> = 0. 

In order to put the theory in a form suitable for the application of the for- 
malism developed in Section 2, we first resort to the Hartree approximation. 
Taking the expectations over the free Gaussian measure with their proper 
haffnian combinatoric factor, the Hartree prescription yields the following 
substitution: 

4 (x,t)^6((0 2 ))^ n 2 (x,t), (36) 

where {{4> 2 ))^J n is the asymptotic, time- independent value of the two-point 
connected Green's function, to be determined self-consistently by equating 
its expression to the equilibrium value {{4> 2 ))^ derived from the dynamics. 

In the Hartree theory the potential V(<f>) is then substituted by an effective 
linear potential 

V cS = ^f<P\ (37) 
where the Hartree effective mass is now given by 

m\ = m 2 + 3A((0 2 »f\ (38) 

The fact that m 2 H is time-independent allows us to solve the eq. fll3|) with 
the substitution m 2 — > m 2 H and to make use of the results given in Subsection 
2.2. The contribution to ((<p 2 {t)))v from the initial configuration is expected 
to be dissipated away within a typical relaxation time r 7 , whereas the con- 
tribution from the stochastic force increases with time reaching for t ^> 1/rj 
the asymptotic value ((0 2 ))y yn - 

From eq. ( |3TD one can easily read off the self-consistency condition as an 
intersection equation: 



«<f» d / n = /(A,«0 2 »r 



rf3 q e-M^ T 



(2tt) 3 \q\ 2 + m 2 H 



rf3 q e -|q| 2 R 2 ? 

(2tt)3 |q| 2 + m2 + 3A((0 2 ))^ n ' 
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An upper bound for ((4> 2 ))v rn then is trivially found observing that 

/(A,((0 2 )) d / n )</(o,((0 2 )) d / n ), (40) 
which, taking into account the results of the previous Section, assures that 

((0 2 )) d / n < 10- 2 mT. (41) 

This inequality is not surprising: being m 2 > 0, a positive value of A has 
the effect to increase the effective mass around = making fluctuations 
more difficult than in the case A = 0. 

Some comments are in order here. First of all, one can expect the Hartree 
approximation to be reliable in the limit of small coupling A, i.e. when 
nonlinear effects can be neglected in first approximation. In this respect, 
the Hartree approximation amounts to consider the scalar field to be the 
modulus of a vector field <p with a large number N of components and self- 
coupling (A/iV) \0\ 4 ]T7 . 



We have considered a nonlinear potential (pq), that does not depend upon 
the volume V of the region where one takes the space average of the variance 
at long times. 

This is strictly correct [[Hj only for regions with typical size R > £ which 



are the ones we are interested in. For R < £, one should modify expression 
( |35D to properly take into account the finite size effects. 

In the next Section we finally address the role of thermal fluctuations 
in the physically interesting case of the electroweak phase transition and 
investigate whether the latter is hindered by the presence and the equilibrium 
dynamics of the subcritical bubbles. 



4. Fluctuations at the electroweak phase transition 
4.1 Small supercooling limit and thin wall approximation 

First order phase transition and critical bubble dynamics in the Standard 
Model have lately been studied in much detail, and it has become increas- 
ingly clear |18|] that for Higgs masses considerably heavier than 60 GeV, the 



electroweak phase transition is only of weakly first order. 
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For a Higgs mass M# > 100 GeV, both perturbative and lattice calcu- 
lations confront technical problems. However, it is conceivable that for such 
large Higgs masses the electroweak phase transition is close to a second order 
and does not proceed by critical bubble formation. 

In this Section we consider a phenomenological Higgs potential for the 
order parameter <ft suitable for a simple description of a first order phase 
transition: 

V(<j>) = \m 2 {T)4? - l -aT^ + ^A0 4 , (42) 

where we have not determined the parameters perturbatively but fit them, 
when needed, according to a recent two-loop determination of the gauge- 
invariant effective potential fl~9|| . 

Most of the dynamical properties of the electroweak phase transition as- 
sociated with the potential Eq. fl4"2|) , such as the smallness of the latent heat, 
the bubble nucleation rate and the size of critical bubbles, have been dis- 
cussed in PDJ. For the purposes of the present paper it suffices to recall only 



some of the results. 

Assuming that there is only little supercooling, as seems to be the case 
for the electroweak phase transition, the bounce action can be written as 

a 2 9 / 2 tt A 3 / 2 

where A(T) = 9Am 2 (T)/(2a 2 T 2 ). The cosmological transition temperature 
is determined from the relation that the Hubble rate equals the transition 
rate oc e~ s / T , yielding S/Tf ~ \n{MpjTj) ~ 150, where Tj is the transition 
temperature. Thus we obtain from Eq. (|43| ) 

a 1 ' 2 

A(7»~ 1 -0.0442^ = 1-5. (44) 

On the other hand, small supercooling implies that 1 — A = S <C 1, i.e. 
a 500 A 3//2 . Solving for A in Eq. ( ^3|) yields the transition temperature Tf. 
One finds 

m 2 (T f ) = ^ A(T/) T 2 . (45) 
The extrema of the potential are given by 

± (T) = |^(l± V / l-8A/9). (46) 
16 



Expanding the potential at the broken minimum <f>+(T) we find 
-6 ee V(<f> + ,T f ) = \m\T f )4>l-^\4>\ = -0.00218^^ + (S 2 ) . (47) 

The height of the barrier is situated at 0_ ~ 0+/2 with V (0_, T c ) = V max = 
a 4 T^/(UA A 3 ), where T c is the temperature at which V(0) = V(<f>+), given 
by the condition m(T c ) 2 = (2 a 2 T 2 /9 A). As T c ~ T/ we may conclude that 
the thin wall approximation is valid if —e/V milx = 0.314 a 1 / 2 /A 3 / 4 <C 1, or 
a 10A 3 / 2 . Thus the small supercooling limit is clearly satisfied if the thin 
wall approximation is valid. 

At T = Tf the system can tunnel from the unbroken phase = to 
the broken phase = 0+ throught the formation of energetically favoured 
critical bubbles with critical radius R c . 

To get the size of the critical bubble we still need the surface tension. 
One easily finds 

*=f W^ea = |^*?- («) 

We define the critical bubble radius by extremizing the bounce action. In 
the thin wall approximation the result is 

A 3 / 4 

if c = i3.4 (49) 
a v i A m\Tf) 

Therefore R c is much larger than the correlation length £(7/) = l/m(Tf) 
at the transition temperature, as it should. It is clear that fluctuations of 
the Higgs scalar field with spatial correlations comparable to the critical 
bubble size may be expected to be important for bubble nucleation. 



4.2 Subcritical bubbles 

Let us first make the general observation that it is the actual transition 
temperature Tf rather than the critical temperature T c which is relevant in 
the study of subcritical bubbles. This is true in the sense that if subcritical 
bubbles are not important at Tf, they will not be at T c either. 

As we shall show, it actually turns out that subcritical bubbles are not 
important even at Tf. This justifies, in retrospect, our choice T = Tf for 
performing the calculations. 



17 



In the case of a weak first order phase transition the critical bubble is 
typically well described by a thin wall approximation, where the configuration 
is by no means gaussian, but has a flat 'highland' (with determined by the 
non-zero minimum of the potential) and a steep slope down to <fi = 0. 

Therefore it seems natural that also a large subcritical bubble should 
resemble the critical one: when its radius R increases, the form of the sub- 
critical bubble should deform smoothly so that, when R = R c , the bubble is 
a critical one. 

The authors of ref. || took into account this observation in their choice 
of the ansatz = 0(0o, R) for the configuration describing the thermal fluc- 
tuations around equilibrium. 

Treating the amplitude 0o (at the core) and the radius R of such fluctu- 
ation as statistical degrees of freedom, they found that the average size (R) 
is much larger than the correlation length £(T/) and that the most probable 
amplitude at the core (0o) 1//2 is much smaller than 0_(T/), the value of the 
field where the potential acquires its maximum. 

It was then concluded M that subcritical bubbles, at least for Higgs 
masses less than about 100 GeV, are not important during the onset of the 
electroweak phase transition. 

This conclusion was drawn also by taking into account the crucial role 
played by thermalization: since the thermalization rate rj for small amplitude 
scalar fluctuations and large spatial size was estimated in the Standard Model 
2lJ to be of order of 10~ 2 T near the critical temperature, i.e. much larger 



than the typical first order transition time, small amplitude fluctuations with 
size larger than about 1/rj were considered to be absent from the mixture of 
subcritical bubbles and not counted in the thermal averages. 

However, we know that the fluctuation-dissipation theorem describes the 
common origin of dissipation and fluctuation: the impact between the sys- 
tem and the surrounding particles of the thermal bath is not limited to a 
systematic friction force responsible for the dissipation of fluctuations, but is 
also characterized by a random force responsible for the fluctuations them- 
selves. This stochastic force is fundamental, as we have seen in Section 2, 
in determining the dynamics of fluctuations with time and expecially in the 
asymptotic regime. 
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From the considerations above, it should be clear that the thermal aver- 
ages (R) and (0q) 1//2 found in ref. || must be respectively regarded as the 
most probable size and amplitude at the core of a subcritical bubble formed 
at t — 0. Let us then imagine that this "most probable" configuration ap- 
pears at t — 0. 

To follow its dynamics, again we linearize the potential applying the 
Hartree approximation: 

3 (x,t) -> 3((0 2 )) d / n 0(x,t), 4 (x,i) -> 6((0 2 >> d /V(x,t). (50) 

The linearized Langevin equation now reads 

0(x, t) - V 2 </> (x, t) + #(x, t) + (x, t) — aT ((0 2 )) J yn = 7 (x, t) , (51) 

where m 2 H is the same in equation (|38|). 

Since we are assuming the electroweak phase transition to be first order, 
which is expected to occur for for light Higgs masses or equivalently A < 1, 
and since thermal fluctuations are expected to be important for a very weak 
first order phase transition, i.e. for a < 1, we expect the Hartree approx- 
imation for the quartic potential fl42[) to be reliable since the coefficients of 
the nonlinear terms in the potential are very small. 

The Hartree approximation does not allow to describe the nonlinear ef- 



fects present in the theory. For instance, Copeland et al. |T6[ have recently 
pointed out that nonlinear scalar field theories can incorporate well localized, 
time-dependent, configurations, called oscillons, which, although unstable, 
can be extremely long-lived, r s ^ ~ (10 3 — 10 4 )m _1 . However, such configura- 
tions can develop resonances only if they form with an initial amplitude 0o a t 
the core above the inflection point 0_, which seems to be not very probable 
M. Moreover, oscillons were studied in absence of any coupling to a thermal 
bath whose presence is expected to reduce considerably their lifetime. 

Equation flSTp is far from being the true Langevin equation describing 
a thermal fluctuation in the hot electroweak theory. Indeed, even in the 
simpler A0 4 theory, the noise is colored and the coupling to the thermal bath 
is multiplicative. 

Also, the viscosity coefficient r\ is not expected to be a constant, but to 
depend upon the typical momentum scale of the thermal fluctuation, i. e. the 
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inverse of its radius, and its amplitude. Up to now the only estimate of rj 
for the Standard Model is that given in ref. and is valid only for small 
amplitude 0o and large size R fluctuations, R > 1 / 77 ~ 10 2 /T. 

Fortunately enough, this is just the range of 0o and R provided in ref. 
[§] for the initial conditions that we have to specify for solving the linearized 
Langevin equation; we therefore expect that assuming a constant n in mo- 
mentum space provides a reliable approximation. 

Again, physical results at long times, or at equilibrium, are obviously 
expected to be independent from the viscosity coefficient 77 and from the 
nature of the thermal noise, even if relaxation times can change significantly. 

Equation (^l|) is easily reduced in a form equivalent to the one holding in 
the free field case by absorbing the constant non-homogenous term through 
the following linear redefinition of the field variable: 

4> (x, t) = <£>(x,i) + K, 



3 \ + {m%/ ((4?))$*)' 



where <^(x, t) solves the Langevin equation (|l"3|) with m 2 — > m 2 H . Note that 
k in the limit of very large ((0 2 ))y yn approaches <p 2 _. 

For definiteness, we fitted our phenomenological potential Eq. (|42|) to the 



two loop result for the effective potential calculated in |T9| for the Higgs mass 
M H = 70 GeV. This yields a ~ 0.048 and A ~ 0.061. One readily verifies 
that we are indeed safely in the thin wall limit. With these parameter values 
(in units of R c ) 

T f = 85.70, 

m 2 (T f ) = 56.78 

r] ~ 10^ 2 T f ~ 0.875 

<j>_(T f ) ~ 24. (53) 

Assuming to be again in the situation rj < 2 m# (we shall verify a posteriori 
that this is indeed the case), we easily read off the solution of the linearized 
Langevin eq. fl5"T|): 



(x,t) = <p(x,t) + K l-e-" t/2 



cos 



?7 2 — 4 m 2 H 
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- e-*A 1 sin V" '- 4 "" t , (54) 
Jrj 2 — A m 2 H 2 y 

where ip(pc,t) does have the same time evolution described in Section 2. 

Dynamics evolves as follows. At £ = a subcritical bubble is formed with 
initial amplitude at the core 0o an d radius R given in ref. ||. For t <C 1/??, 
{{4> 2 (t)))v none of the modes of the initial configuration have yet relaxed 
to zero, even if they start to decrease as £ 2 . In the meanwhile, besides the 
noise term ((<p 2 {t)))y which increases as £ 3 , a new term begins to contribute to 
(((j) 2 ))v, i-e. the term proportional to k 2 increasing as £ 2 . For£ ^> 1/rj, all the 
modes have relaxed towards their equilibrium value. The initial configuration 
has disappeared, and ((4> 2 ))'y n is the sum of the noise term contribution and 
the K 2 -term contribution, so that the self-consistency equation for ((4 >2 ))'v n 
reads 



«0 2 »r = ;(«0 2 »r) 



-2. / i!SU -lql 2 * 2 T 



K + ' |q| 2 + m2 + 3A((0 2 )) d / n ' (55) 

In Fig. 1 we present the behaviour of the function / (((0 2 ))y yn ) (dashed line) 
in units of 1/R C for A = 0.061, a = 0.048 and R = £(Tf), i.e. the variance 
is calculated in a typical correlation volume. Note that in the limit of very 
large ((0 2 ))^ yn , / tends toward (j> 2 _(T f ) ~ (24/R c ) 2 . 

Fig. 2 gives the same function for the two cases R = £(Tf) and R = 
2 £(Ty). The points where the dashed lines intersect the solid line, repre- 
senting the left hand side of eq. (|55|) , denote the values of ((<p 2 ))y' n which 
satisfy the self-consistency equation in the two different cases: when testing 
the fluctuation in a correaltion volume = (47r/3)£ 3 , we have 



* (2.7/R c ) « (j>-{Tf) ~ (24/i? c ); (56) 



if we test the fluctuation in a volume larger than V^, the variance is even 
smaller, 



'((^)) d /^(l. 2 / J R c ) (57) 
for R = 2£(Tf). Note also that our starting assumption r] < 2m H is verified. 
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One could naively think that the extra energy due to a nonvanishing vari- 
ance would facilitate barrier penetration for the formation of critical bubbles 
or possibly even invalidate our thin wall approach. 

We can verify that this is not the case by considering the fluctuation 
energy about = 0. For R = £(T"/) it is given approximatly by (in units of 

E a ~ ^m 2 (T / )((0 2 ))^ yn ~ 211 « -e ~ 5 x 10 3 , (58) 

where we have made use of the fact that + (Tj) ~ (48/i2 c ). Thus, in terms of 
the barrier penetration, thermal fluctuations represent only a minor correc- 
tion and they do not have remarkable effect on the thin wall approximation. 

From these considerations and from the fact that the variance of the 
thermal fluctuations \J ((0 2 ))y yn is clearly smaller than the inflection point 
0_(T/), we can conclude that subcritical bubbles do not affect the nucleation 
of critical bubbles in an appreciably way. 



5. Conclusions and perspectives for future work 

In this paper we have investigated the role played by subcritical bubbles 
at the onset of the electroweak phase transition. 

Treating the the configuration (x, t) modelling such thermal fluctuation 
as a stochastic variable, we have described its dynamics by a phenomenologi- 
cal Langevin equation for different models. This approach allows to properly 
take into account both effects of the thermal bath on the system: a sys- 
tematic dissipative force, which tends to erase out any initial subcritical 
configuration, and a random stochastic force responsible for fluctuations. 

Following the evolution with time of the variance {{<p 2 {t)))v i n a given 
volume V, we have shown that the contribution to it from any initial sub- 
critical configuration is quickly damped away and, in the limit of long times, 
{{4> 2 (t)))v approaches its equilibrium value provided by the stochastic force 
and independent from the viscosity coefficient r], in agreement with the 
fluctuation-dissipation theorem. 

In the most interesting case of thermal fluctuations during the electroweak 
phase transition, we have made use of a self-consistent Hartree approximation 
expected to give good results in the weak coupling limit which we believe to 
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hold for a weak first order phase transition. A more correct way of proceeding 
might have been to make use of the variance ((<p 2 {t)))v i n the substitution 
( |50D instead of ((0 2 )}y yn - However, since the system reaches very quickly the 
equilibrium, we are confident that our approximation is reliable. 

We have shown that thermal fluctuations do not affect the nucleation 
of critical bubbles at the tunneling temperature Tf and from this fact we 
conclude that electroweak baryogenesis scenarios associated to a weak first 
order electroweak phase transition remain a viable possibility to explain the 
primordial baryon asymmetry in the Universe. 

Very recently Bettencourt |J has computed the probability at equilibrium 
for the fluctuations of the Standard Model Higgs field, averaged over a corre- 
lation volume, to exceed the inflection point 0_ at T ~ Tf and for a physical 
Higgs mass M# = 70 GeV, the same value we have used in our paper. 

He estimated such a probability to be very small, around 1.27%, and to 
decrease when the testing volume V is increased. These estimates and be- 
haviour with V are in complete agreement with our results for the asymptotic 
value of ((<f) 2 (t)))v whose evolution with time is also given here. 

The phenomenological Langevin equation used in our paper is not the 
true equation describing a thermal fluctuation in the hot electroweak theory. 

In this respect, our work should be regarded as a first step towards a com- 
plete description of the dynamics of thermal fluctuations at the electroweak 
phase transition by a more complicated Langevin equation which properly 
describes the coupling of the stochastic field <fi (x, t) to the other degrees of 
freedom. 

Nevertheless, we expect our physical results to be reliable at long times 
when the details on the nature of the thermal noise and the viscosity co- 
efficient 1] become unessential. On the other hand, relaxation time-scales 
crucially depend on the nature of the stochastic force and the strength of 
dissipation and their complete knowledge is needed to decide if degrees of 
freedom other than (x, t) are in equilibrium or not inside the subcritical 
bubble. 

One may realize that this is a crucial question by reminding that the 
effective potential (|42"l) , used to describe the free energy associated to the 
fluctuations, is usually obtained at 1-loop integrating out all the degrees of 
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freedom of the theory other than (f) (x, t), i.e. fermions and gauge bosons. 

In performing such a calculation, it is assumed that fermions and gauge 
bosons do have equilibrium distributions with a (x, t) background depen- 
dent mass. This is true only if their interaction times with the background 
<fi (x, t) are much smaller than the typical lifetime of the subcritical bubble. 

Since this condition seems not be satisfied, a full non-equilibrium ap- 
proach is needed. Such an approach might lead to unexpected results as far 
as the issue of thermal fluctuations during the electroweak phase transition 
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Figure Captions 

Figure 1: The plot of the function / (((0 2 ))y yn ) (dashed line) in units 
of l/R 2 c for A = 0.061, a = 0.048 and R = £(7)). The solid line represents 
the inflection point (f) 2 _(Tf). 

Figure 1: The plot of the function / (((0 2 ))y yn ) (dashed lines) in units 
of \jR\ for the same values of A and a given in Figure 1 and for R = £(T/) 
and R = 2£(T/) respectively. The intersection points between the dashed 
lines and the solid line represent the values of ((0 2 ))y yn satisfying the self- 
consistency equation in the two cases, see the text. 
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